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Abstract: We derive closed expressions for the mean and variance of the (mini-)jet 
multiplicity distribution in hard scattering processes at low x. Here (mini-)jets are 
defined as those due to initial-state radiation of gluons with transverse momenta 
greater than some resolution scale /i R , where A 2 <C /i| <C Q 2 , A being the intrinsic 
QCD scale and Q the momentum transfer scale of the hard scattering. Our results 
are valid to leading order in log(l/x) but include all sub-leading logarithms of Q 2 / /x R . 
As an illustration, we predict the mini-jet multiplicity in Higgs boson production at 
the Large Hadron Collider. 

Keywords: Deep Inelastic Scattering, QCD, Jets, Hadronic Colliders. 



*Work supported in part by the UK Particle Physics and Astronomy Research Council and by 
the EU Fourth Framework Programme 'Training and Mobility of Researchers', Network 'Quantum 
Chromodynamics and the Deep Structure of Elementary Particles', contract FMRX-CT98-0194 
(DG 12 - MIHT). 



Contents 

1. Introduction 1 

2. Jet-rate generating function 2 

3. Jet multiplicity moments 5 

4. Discussion of the jet multiplicity moments 7 

5. Mini-jet multiplicity in Higgs production at LHC 10 

6. Conclusions 12 
A. Lipatov anomalous dimension and characteristic function 12 



1. Introduction 

The prediction of jet multiplicities in hard scattering processes represents an impor- 
tant, well-defined challenge to perturbative QCD. Especially challenging is the calcu- 
lation of the jet multiplicity at small x, that is, in processes where the hard-scattering 
scale, Q 2 , is much less than the overall collision energy-squared s (x ~ Q 2 /s <C 1). 
Such processes include not only the familiar case of deep inelastic lepton scattering 
at low values of the Bjorken variable, but also Higgs and gauge boson production 
at high-energy hadron colliders (see e.g. ref. [1]). A large associated jet multiplicity 
could pose serious triggering and background problems for experiments searching for 
new physics at the future Large Hadron Collider (LHC). 

The particular problem at small x is that higher-order perturbative contributions 
may become enhanced by factors of ln(l/x). As a first approximation, one can 
consider only those terms that are enhanced by a factor of m(l/a;) for each power 
of the strong coupling a s . In this approximation, the fully inclusive deep inelastic 
structure functions satisfy the leading-order Balitskii-Fadin-Kuraev-Lipatov (BFKL) 
equation [2, 3]. 

When computing the associated multi-jet production rates at small x [4, 5, 6], 
one encounters another potentially large logarithmic factor, namely ha.(Q 2 //j%) where 
/i R is the minimum transverse momentum that a jet must have in order to be resolved. 
For mini-jets, by which we mean jets with transverse momenta that are relatively 
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low compared with the hard scattering scale, it clearly becomes important to resum 
terms enhanced by factors of ln(Q 2 //i|) as well as ln(l/x). 

In ref. [5], such a resummation was performed in the double-logarithmic (DL) ap- 
proximation, i.e. keeping only terms of the form [a s ln(l/x) ln(Q 2 //i|)] n to all orders. 
The result was remarkably simple, namely that the mean and variance (mean square 
fluctuation) of the jet multiplicity are respectively quadratic and cubic functions of 
ln(Q 2 //i|), with simple coefficients. 

In ref. [6] we extended the results of ref. [5] concerning fixed-multiplicity jet rates 
to single logarithmic (SL) accuracy, i.e. keeping all terms [a s ln(l/:r)] n [\n(Q 2 / fi^)] m 
with < m < n. We also showed that an iterative procedure could be used to obtain 
the jet multiplicity moments to any desired order in a s , to SL accuracy. However, 
we were not able to derive closed expressions that resum the enhanced terms in the 
multiplicity moments to all orders. That is done in the present paper. 

An alternative approach to hard scattering at small x is the CCFM equation 
[7, 8] based on angular ordering of gluon emissions. Remarkably, it can be shown 
[9] that for sufficiently inclusive observables the CCFM formalism leads to the same 
results as the BFKL equation. The jet multiplicities calculated in the present paper 
fall into this class of observables. The BFKL formalism turns out to be technically 
simpler and therefore we adopt it here. 

The paper is organized as follows. First, in sect. 2, we solve the equation derived 
in ref. [6] for the jet-rate generating function, by expressing the solution as an inverse 
Laplace transform. In sect. 3 we obtain the jet multiplicity moments by saddle-point 
evaluation of the inverse Laplace transform. The results are discussed in sect. 4. We 
find that the simple quadratic and cubic forms of the mean and variance found in the 
DL approximation remain valid to SL precision. The coefficients, while not particu- 
larly simple, can be expressed straightforwardly in terms of the Lipatov anomalous 
dimension [2] and its derivatives. 

As an illustrative application of our results, in sect. 5 we calculate the mean 
and variance of the mini-jet multiplicity in Higgs boson production at the LHC, as 
functions of the Higgs mass. 

2. Jet-rate generating function 

We write the contribution to the gluon structure function at scale Q 2 , F(x,Q 2 ), in 
which r final-state mini-jets are resolved with transverse momenta greater than /i R , 
in the form 

F (rjct) (a;, Q 2 , /4) = F(x, /4) <8> G ( - r \x, T) = (' -F(z, &)& r \x/z,T) , (2.1) 

Jx Z 

where the coefficient is a function of x and T = ln(Q 2 //-4) to be determined . 
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The Mellin transformation 



/«(•••) = / dx aff(x,...) , 
Jo 



with inverse 



lJ cd . x -^u...), 



(2.2) 



(2.3) 



where the contour C is parallel to the imaginary axis and to the right of all singu- 
larities of the integrand, converts the convolution in eq. (2.1) into a simple product: 

F^ ct \Q 2 ^l) = FM)G^\T). (2.4) 

Furthermore the evolution of the structure function in u;-space is simply given by 



F w iQ 2 ) = exp[ 7L (a s /cu)T]i^(/4) , 
where 7 L is the Lipatov anomalous dimension: 



7L (a s /cu) = ^ + 2C(3) 



UJ 



UJ 



+ 2C(5) 



0: t 



UJ 



+ 12[C(3)f 



UJ 



+ 



which is the solution of the equation 

uj = -a s [27 E + ^(7) + ipil - 7)] = a s xil) , 



(2.5) 



(2.6) 



(2.7) 



where a s = 3a s /n, ip being the digamma function and 7 E = — tp(l) the Euler con- 
stant. 

Introducing the jet-rate generating function 



GUu,T) = Y,u r G ( r\T), 



(2.8) 



r=0 



we can write the r-jet rate, i. e. the fraction of events with r mini-jets, as 

1 d r 



r! du r 



/?„(«, T) 



u=0 



where 



Kiu,T) = exp[- 7L (a s /cu)T] G W («,T) . 
We showed in ref. [6] that G w (w,T) is given by 



(2.9) 



(2.10) 



G w (u,T) = l+ / dtg u iu,t,T), (2.11) 
where the unintegrated function g u iu,t,T) satisfies the integro-differential equation 



0« 



1 + / dt' (u, max{t, t'}, T) + 2 V C(2m + 1) 



(2.12) 
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Rearranging terms and differentiating with respect to t, this gives 



d_ 

dt 



1 - (u-l)t 



\ a 

guj > = -u— 

J 



s 
CO 



i-2£c(2m + i; 



m=l 



Q2m+1 

dt 2m+1 



(2.13) 



To find g w we first recall the solution in double-logarithmic approximation [5], 



g^( U ,t,T)=U-^ 



l-(u-l)*— 



1 - (u - 1)T— 



(2.14) 



Notice that this expression vanishes, together with all its t-derivatives, at the point 
t — to — u/[a s (u — 1)] as u — > 1 + . We shall assume this is also true to SL accuracy, 
for some finite value of t , and then verify that the solution is consistent with per- 
turbation theory. The equation (2.13) can now be solved by introducing a Laplace 
transformation of the form 



t 



dtg u (u,t,T) e 7 * , 



with inverse 



9u(u,t,T) = ^-r J^d'jg U! (u,-f,T) e 



■it 



(2.15) 



(2.16) 



where the contour V is parallel to the imaginary axis and to the right of all singular- 
ities of the integrand. This gives 



^(m,7,T) = ^(w,7o,T)e 



(2.17) 



where 



u 



u 



1 Jj 



LO 



a s u 



- x(Y) 



u 



u-l 



\a s u 



+ 27 E ) (7 - 70) - In 



r(i- 7 )r( 7 o) 
r( 7 )r(i- 7o ) 



(2.18) 



X being the Lipatov characteristic function, see eq. (2.7). To determine gu(u, 70, T) 
we make use of the boundary condition provided by eq. (2.12), 



g u (u,T,T) = 

which tells us that 
1 



to + a s T 



l + gJ + 2 £ C(2m+1) 



m=l 



d 2m g u 

Q t 2m 



Jt=T 



2m 



-gu,(uno,T) / c?7 



00 «-L 1 / x 

T+-- X (l) 

a s u u 7 



e -7T+<Mu, 7 ) — 1 _ 



(2.19) 



(2.20) 



Deleting a total derivative from the integrand, since endpoint contributions must 
vanish, we obtain simply 



(2.21) 
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where 



Hence 



l„(u 



■ T > = /r 



c?7 



-7T+<M«:7) 



r 7 



so that finally 



g w (u,t,T) = / / d-fi 
I u (u,T) Jr 

GUu,T) = 1 + [ T dtgUu,t,T) = ^4 

JO 1uj{U,1) 



(2.22) 
(2.23) 

(2.24) 



3. Jet multiplicity moments 



The moments of the jet multiplicity distribution are obtained by successive differen- 
tiation of the generating function at u — 1: 



r(r — 1) . . . (r — s + 1) 



<9w s 



exp[-7 L (a s /^)T] 



d s G u 



u=l 



du s 



(3.1) 



u=l 



To compute these quantities to all orders from the expression (2.24) we can use 
the method of steepest descent. It is clear from eq. (2.18) that the saddle point of 
<f>u>( u ,l) occurs at 7 = 7 S where u = a s -ux(7s), or in other words, from eq. (2.7), 
7s — 7l(«s m / luJ ) ■ Choosing for convenience 7 = 7 S , we have 



&>(u,7) 



I E^(7-7 s rX { "- 1) (7 s ) 



u 



(3.2) 



n=2 



Now we have to evaluate integrals of the form 1 

oo j{pi) 



f d 7 /(7)e*- (t, ' 7) = / rf 7 E ^V(7-7s) 
Jr Jr ^ Q ml 



= {fs + 



(u-1) A „ 1 X " / 



_ fll _ a.s ri 

9 2 x 



+ 



u-1 



2 r 



ill 



5 y" 1y"' 5v" 

^_ f/l/l ^_ As fill _ As fll _|_ _ As /•// 

8 /s 12 xi Js 4 X '/ s+ 8x'r s 



8 X ' Js 



3 X? L 8 Vx' s 



+ . 



7) 



(3.3) 



where f s = /(7 S ), x s = x'(7s), etc., and the dots represent terms of higher order in 
(u — 1). To express everything in terms of /(7 L ), x'(1l), etc., where 7 L = 7 l (q; s /c<j), 
we use the results in the Appendix to write 

7s = 7 L +7a--iA + k>-i) 2 



= 7l - (« - 1)^ + (« - I) 2 (- - ^C^l + 



/3 



(3.4) 



1 The saddle-point method has been applied to similar integrals in ref. [10]. 
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where now \ — x(7l), etc., and we can use this to expand eq. (3.3). 
In the present case we have 

/(7) = -e~" T ■ (3.5) 

7 

Hence from eqs. (2.24), (3.1) and (3.3) we obtain for the mean number of jets, 



1 



r u , = 



X 



X' \7l 2 X 
and for the variance a 2 , = r 2 w — r 2 ,, 



1 + X T 



2x r 



rj~i2 



(3.6) 



cr„, = — 




x 2 x" 



X 



12 



3X" 

ilx' 2 



2 X " 2 


5%" 3 


7lX' 3 


4x' 4 


X"" ' 


T 


x" 2 


x'" 


X' 3 


2X' 2 



rp2 



(3.7) 



Eqs. (3.6) and (3.7) are convenient for numerical calculation, since they only 
involve the evaluation of polygamma functions. To compare with perturbative ex- 
pansions, we can use the results in the Appendix to rewrite them in terms of 7 L and 
its derivatives: 



o,. 



u \7l 2 7 ^ 



a Sr ^ 1 , / a s 



T 



'_2 + 

v 7l 7l L 2 7 ^ 



7f 7f 2tT 



7T 

7[, 



+ 



-1 + + — 
2 



1 



+ 3 ^ + ^ 



^->3 



(3.J 




a; 



(3.9) 
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Using the perturbative expansion (2.6) for 7 L , we find 

uj 2 V uj J \uj J 



+4C(3) (^)^ 2 - 8C(5) (^) 6 T + • • • (3.10) 



and 



ol = + V T 2 + | f ^ 3 T 3 - 2 C (3) f ^ 4 T 



+12C(3)(^) 5 T 2 -(8C(5)-|C(3)T 2 )(^)°T + ... , (3.11) 
which do indeed agree with the expansions obtained in ref. [6]. 

4. Discussion of the jet multiplicity moments 

We see from eqs. (3.6) and (3.7) that it remains true to all orders to SL precision that, 
as in the DL approximation [5], the mean number of jets is a quadratic function of T 
while the variance is a cubic function of T. Thus the distribution of jet multiplicity 
at small x and large T remains narrow, in the sense that its r.m.s. width increases 
less rapidly than its mean with increasing T. 

It is interesting to see how rapidly the perturbative expansions for the mean 
number of mini-jets and its variance converge. We first introduce the notation 

^ 0l (f) r + 02 (f)V (,!) 



and 

^(f) r + ^) V + <*(f) V <"> 

Fig. 1 shows the coefficients a; (i = 1,2) as functions of a s /uj, including the exact 
dependence according to eq. (3.6) as well as the perturbative expressions (3.8) ex- 
panded to different orders. The coefficients bi (i = 1,2,3) are shown in fig. 2. 



The coefficients in the mean jet multiplicity f u are extremely well approximated 
by the perturbative expressions for small values of a s /uj, up to 0.3. In this region, the 
coefficients remain close to the (constant) values predicted by the double-logarithmic 
approximation. Above a s /uj ~ 0.3, where the coefficients start to diverge, the per- 
turbative expansion converges rather slowly. A similar behaviour is observed for the 
coefficients in the variance <r 2 . The perturbative expansion converges very rapidly 
for small a s /uo, in this case up to ~ 0.25. Above a s /uo ~ 0.3, the convergence is 
again rather poor. 
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Figure 1: Coefficients in the mean jet multiplicity r w as functions of a s /to, as defined in 
(4.1). The solid line represents the exact coefficients as calculated from (3.6), the other 
lines are obtained by expanding (3.8) to 10th, 14th, and 15th order in a s /u. 



The singularities in the coefficients arise from the fact that eq. (2.7) does not 
have real solutions for a s /uo > (41n2) _1 = 0.361. This bound will be raised [10] 
when NLO corrections to the BFKL equation [11, 12] are included. We therefore 
expect higher order corrections to the BFKL equation to have a significant effect on 
the jet multiplicity moments for values of a s /oo above ~ 0.3. 

For practical applications, we need to compute the multiplicity moments as func- 
tions of x. This is done most conveniently using the perturbative expansions (3.10) 
and (3.11), extended to sufficiently high order, and performing the inverse Mellin 
transformation term by term using 



1 

27ri Jc 



du x 



-LO-l 



CO 



a s [q s ln(l/a;)] 
x (n — 1)! 



n— 1 



(4.3) 



The factorial suppression of high orders in x-space removes the convergence problems 
that were encountered in cu-space. 

To illustrate the behaviour in rr-space, we consider the mini-jet multiplicity as- 
sociated with pointlike scattering on the gluonic component of the proton at small 
x. The mean number of mini-jets considered as a function of x is given by 



n(x) = 



F(x,Q 2 ) (g)r(x) 
F(x,Q>) 



(4.4) 



where r(x) is the inverse Mellin transform of r u , obtained by means of eq. (4.3). 
Similarly the dispersion a n in the mean number of mini-jets as a function of x is 



2 F(x,Q 2 )®r 2 (x) 2 
= F(x,Q 2 ) HX)] 



(4.5) 
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To perform the convolution in eqs. (4.4), (4.5) we used the leading-order MRST 
gluon distribution [13] . Fig. 3 shows the resulting x-dependence of the mean number 
of mini-jets n(x) and its dispersion cr n (x) for Q = 100 GeV and two different values 
of the resolution scale fi R . We see that the results remain close to the predictions 
of the double-logarithmic approximation at all but the very smallest values of x and 
fi R . The convergence of the perturbation series in x-space is so good that there is no 
visible difference between the results at 5th and 15th order for the ranges of x and 
/i R shown. 

We should emphasise that the results in fig. 3 are not directly comparable with 
data on deep inelastic lepton scattering (DIS). To make quantitative predictions for 
DIS one would need to perform a further convolution with an impact factor that 
represents the "quark box" coupling of the gluon to the virtual photon. One should 
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Figure 3: The mean and dispersion of the mini-jet multiplicity as functions of x at 
Q = 100 GeV for two different resolution scales /j, r . The solid lines show the SL results 
including the 15th order in perturbation theory, the dashed lines correspond to the DL 
approximation. 



also take into account the production of mini-jets from the quark box, together with 
non-perturbative effects such as mini-jet emission from the proton remnant. 



5. Mini-jet multiplicity in Higgs production at LHC 

The dominant production process for a standard model Higgs boson at the LHC 
is expected to be gluon-gluon fusion. This process is illustrated in fig. 4. The 



Pi 
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Figure 4: Higgs boson production by gluon-gluon fusion. 



production cross section for a Higgs boson of mass M and rapidity y by gluon-gluon 
fusion in proton-proton collisions at centre-of mass energy y/s takes the form [1] 

^ = F(xi, M 2 ) F(x 2 , M 2 ) C{M 2 ) , (5.1) 
ay 
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(5.2) 



where 

M y M 

%i = —=e y , x 2 = —= e y , 
V s V s 

and for LHC \fs = 14 TeV. C represents the gg — > H vertex, which is perturbatively 
calculable as an intermediate top-quark loop. In the mean number of mini-jets and 
its dispersion, however, C cancels and we do not need its detailed form. Thus we can 
compute the associated mini-jet multiplicity in Higgs production quite simply using 
the machinery developed in the previous section. 

We will consider central production of the Higgs boson (y = 0), and therefore 
we have x\ = x 2 = x where x = M / y/s. Since the gluon emissions in the upper and 
lower parts of fig. 4 are independent, we can simply add the numbers n\ = n(xi) 
and n 2 = n(x 2 ) of mini-jets, and the mean multiplicity N of associated mini-jets 
becomes 2 

N(x) =m + n 2 = 2n{x) , (5.3) 

where n(x) can be calculated as in (4.4) after replacing Q 2 by M 2 . Similarly, the 
variance is 



a 2 N (x)=<j 2 n (x 1 )+<7 2 n (x 2 )=2<j 2 n (x) . 



(5.4) 



Again, a 2 can be obtained as in eq. (4.5). 

We have calculated the dependence of and a N on the Higgs mass M and our 
numerical results are shown in fig. 5. Here, the DL results give an excellent approx- 
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Figure 5: The mean value and dispersion of the number of (mini-)jets in central Higgs 
production at LHC for two different resolution scales fi R . Solid lines show the SL results up 
to the 15th order in perturbation theory, dashed lines correspond to the DL approximation. 



imation and the SL terms are less significant. We see that the mini-jet multiplicity 
and its dispersion are rather insensitive to the Higgs mass at the energy of the LHC. 



2 Again we do not count any jets emerging from the proton remnants. 
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6. Conclusions 



We have derived expressions for the mean and variance of the jet multiplicity distri- 
bution at small x, including resummation of leading logarithms of x and all (leading 
and sub- leading) logarithms of Q 2 / Our results have been derived using the BFKL 
equation, but are expected to hold in the CCFM formalism as well. 

Considered as functions of uj, the moment variable Mellin-conjugate to x, our 
expressions exhibit bad behaviour at large values of a s /uj, which is associated with the 
singularity of the leading-order Lipatov anomalous dimension 7 L at a s /uj = (4 In 2) _1 . 
We would expect this behaviour to be modified strongly by higher-order corrections. 
Although the next-to-leading corrections to 7 L are known, a full calculation of the 
corresponding corrections to the associated jet multiplicity has not been performed 
and would appear much more difficult. 

In x-space we have only been able to find the jet multiplicity and its variance as 
perturbative expressions. These are obtained after expanding the closed expressions 
in cj-space to sufficiently high order and applying the inverse Mellin transformation 
term by term. For all practical applications this is sufficient since the perturbative 
series in x-space turns out to be very rapidly convergent. 

The multiplicity of mini-jets at small x is an observable that can also be computed 
in Monte Carlo simulations of BFKL dynamics [14]. It would be interesting to 
compare the corresponding results with our analytic results. 

A. Lipatov anomalous dimension and characteristic function 



For changing from 7 L and its derivatives to x and its derivatives we need the following 
relations, which are easily obtained by repeated differentiation of the first: 

UJ 

a s 

' uj 

1 

^s/ 7l 



x(7l) 
X'(7 L ) 
X"(7 L ) 
X'"(7 L ) 
X W '(7 L ) 



2 1 



UJ \ " 1 / ^ / C^S ;/ 



Vtts/ 7l 

f-V- 

U s 7 7[ 7 



-67l - 37 L - 6— 7l7l + I — 7l7 

UJ \ CO 



(A.l) 
(A.2) 
(A.3) 
(A.4) 



24 7 f + 36%^' + 30 p 

UJ \UJ 



« 2 - 



«S ^ ,2 /// 

uj ■ 7l 7l 



3 



a. 



+15,^ Tf-lO^J TM'+^I ill 



a. 



3 



am, 

L IL 



Conversely, the derivatives of 7 L are given in terms of x(7l) by 

x 



.(A.5) 



(A.6) 
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< = §3 (ixx' 2 - x") 

<' = 4 [" 6 xV + 6 xx'V - 3x" 2 + x'x ! 



(A.7) 



(A.8) 



7 r = \ [24 X \' 6 - 15 X " 3 + 36 XX 'Y' 2 - 36 X YY - 12xx'\' 



+Wx"x" - x' 2 x' 



(A.9) 
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